Effective squeezing enhancement via measurement-induced non-Gaussian operation 

and its application to dense coding scheme 
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We study the measurement-induced non-Gaussian operation on the single- and two-mode Gaus- 
sian squeezed vacuum states with beam splitters and on-off type photon detectors, with which 
mixed non-Gaussian states are generally obtained in the conditional process. It is known that the 
entanglement can be enhanced via this non-Gaussian operation on the two-mode squeezed vacuum 
state. We show that, in the range of practical squeezing parameters, the conditional outputs are 
still close to Gaussian states, but their second order variances of quantum fluctuations and correla- 
tions are effectively suppressed and enhanced, respectively. To investigate an operational meaning 
of these states, especially entangled states, we also evaluate the quantum dense coding scheme from 
the viewpoint of the mutual information, and we show that non-Gaussian entangled state can be 
advantageous compared with the original two-mode squeezed state. 

PACS numbers: 03.67.Hk, 03.67.Mn, 42.50.Dv 
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I. INTRODUCTION 

Non-classical optical Gaussian states, such as single- 
or two-mode squeezed vacuum states, play essential roles 
in continuous variable (CV) quantum information tech- 
nology. These states have already becn implcmcnted in 
labs and manipulating these states with linear operations 
including beam splitters, phase shifters, displacements, 
and homodyne measurements, various quantum informa- 
tion protocols have been demonstrated: teleportation pj , 
dense coding 0] , and entanglement swapping Q . Mathe- 
matically, the operations by these linear òptics tools and 
arbitrary squeezing operations are described by linear 
and bilincar Hamiltonians and classified as Gaussian op- 
erations which transform Gaussian states into Gaussian 
states. 

A class of Gaussian operation is, however, obviously a 
part of the class of universal quantum operations. Re- 
cent theoretical investigation has shown the limitation of 
Gaussian operations in quantum information processing, 
e.g. entanglement distillation of Gaussian input to Gaus- 
sian output is impossible |j. LJ. |i| anc h more generally, 
quantum information protocols consisting of only Gaus- 
sian operations can be simulated classically 0. There- 
fore, implementation of non-Gaussian operation, that is 
the operations accessible to the outside of the Gaussian 
domain, would be crucial to extract ultimate potential 
of quantum information theory. In addition, a delight- 
ful theoretical result is that, in principle, arbitrary op- 
erations can be implemented by combining one of non- 
Gaussian operations with suitable Gaussian operations 
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Although, at present, even the cúbic nonlinearity is 
hard to realize on the level of single photon, there is an 
alternative idea, called the measurement-induced nonlin- 
earity, where effective nonlinearity is associated with non- 
Gaussian measurement such as photon counting A 
simplest example of such operations has been theoreti- 
cally investigated, where photons in non-classical Gaus- 
sian states are subtracted by low reflectance beam split- 
ters and photon counters. It was proposed that one can 
generate Schròdinger cat-like state from a single-mode 
squeezed vacuum by subtracting photons ^(j- Recently, 
nonclassicality of this cat-like state has been investigated 
with respect to negativity of Wigner function, in consid- 
eration of some experimental parameters 

Moreover, it is predicted that the photon subtraction 
from two-mode squeezed vacuum can increase entangle- 
ment [Ï2l IÏ3l IÏ4||. In ideal situation, i.e. perfect photon 
number counting and lossless setup, the output is always 
pure and one can uniquely quantify the increase of entan- 
glement by the von Neumann entropy of a partial system. 
In practical situation, however, imperfections should be 
taken into account. The main one is the imperfection 
of photon detector. It is still difficult to distinguish the 
photon number precisely. Currently available type of de- 
tector is the on-off type detector (e.g. avalanche pho- 
todiodes in Geiger mode operation) which discriminatc 
only between the vacuum and the presence of photons 
with finite quantum efficiency and nonzero dark counts. 
This type of detector suffices for some purposes. In fact, 
the first observation of non-Gaussian statistics due to the 
photon subtraction from a single-mode squeezed vacuum 
was demonstrated using such a type of device ^íj| . 

A serious restriction due to the on-off type resolution is 
that the detector projeets the original state into a mixed 
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state. For instance, in the case where one of the two 
mode entangled state is measured by the on-off detector, 
the state of the other mode is projected into a mixed 
state. Similarly, the photon subtraction with on-off de- 
tector reduces the pure two-mode squeezed state to the 
mixed state. One cannot, therefore, apply the von Neu- 
mann entropy to quantify their entanglement. In this di- 
rection, some operational measures have been exploited 
theoretically to characterize the non-Gaussian mixed en- 
tangled state. They are, for example, the improvement 
of the teleportation fidelity 0] and the non-locality due 
to the violation of Bell type inequality [Ï7L lïl llH l20| . 

In this paper, we consider the photon subtraction 
scheme consisting of two on-off detectors. The scheme 
is applied to both single- and two-mode squeezed vacua 
including realistic parameters of possible imperfections. 
Special attention is paid to that the density operators 
conditioned by the on-off detector can be represented in 
terms of the sum of the three kinds of Gaussian states. 
Therefore photon subtracted non-Gaussian states may 
still include Gaussian nature more or less. 

Based on this, we address the following two points. 
First, we discuss the validity of the second order variance 
as a measure of the performance of the photon subtracted 
state. Experimentally, the evaluation of the states by the 
second order variance of quantum fluctuations or corre- 
lations are easier and more aceurate than the evaluation 
based on full reconstruction of the states such as quan- 
tum tomography. When two on-off detectors are applied 
to a single-mode squeezed vacuum, one can condition- 
ally generate a plus-cat-like state which is often squeezed 
than the input. It is shown that in realistic squeezing 
regime, the squeezing of the photon subtracted state is 
higher than that of the input. We also show that photon 
subtractions from two-mode squeezed vacuum greatly en- 
hance the second order quantum correlation of the output 
in practical parameter regime. 

Second, we apply the photon subtracted entangled 
state to quantum dense coding |22|. This is an 
entanglement-assisted coding to transmit classical infor- 
mation which can attain a larger capacity than that with- 
out entanglement |2jj. Therefore the increase of the ca- 
pacity, more precisely, the mutual information for a spec- 
ified measurement (e.g. the Bell measurement), can be 
an alternative operational measure for the photon sub- 
tracted entangled state. The mutual information is cal- 
culated by determining the channel matrix between the 
input signals and the measurement outeomes. This mea- 
sure can be a stringent figure of mèrit for the system in 
the sense that the gain usually vanishes even with rela- 
tively small imperfections. This also specifies the asymp- 
totic rate of transmission when the múltiple use of the 
channels is considered. So it would be worth considering 
such kind of information theoretic measure to quantify 
the non-Gaussian state. 

This paper is organized as follows: In Sec. [n] we dis- 
cuss the measurement-induced non-Gaussian operation 
on single-mode squeezed vacuum state, and its equiva- 



lence to Schròdinger cat-like state generation in Ref. ^(| . 
In Sec. Ulli we discuss the two-mode case. The gener- 
ated non-Gaussian entangled state is applied to the dense 
coding scheme in Sec. EO In the Secs. HEU and lVlTl 
we give analyses with consideration of practical parame- 
ters of imperfections. The last section IVllll is devoted to 
discussion and conclusion. 



II. NON-GAUSSIAN OPERATION ON 
SINGLE-MODE SQUEEZED VACUUM STATE 

We ffrst consider the non-Gaussian operation on the 
single-mode squeezed vacuum state in the ideal situation. 
The schematic is shown in Fig. ^ The target mode is 
denoted by path A. For the later extension to the two 
mode case, we consider another mode B whose initial 
state is the vacuum state. The initial state of mode A is 
the squeezed vacuum state 



\r) k = S k (r)\0), 
where S k {r) is the squeezing operator of mode k, 
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and r is the squeezing parameter, proportional to the 
second-order susceptibility and the thickness of the non- 
linear optical crystal and the pump intensity. For this 
state, the uncertainty is reduced in terms of the x— 
quadrature corresponding to ip — in the following 
quadrature operator 
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The input squeezed state is divided into two modes A 
and B with a beam splitter of transmittance t. Mode C 
(D) is then tapped from mode A (B) with a beam splitter 
of high transmittance T. The resulting four-mode state 
is then 
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is the beam splitting operator, and the parameters 4> and 
9 are related with the transmittanecs t and T as 
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respectively. We consider a balanced interferometer (r — 
0.5). Modes C and D are led to the on-off type pho- 
ton detector. The probability operator valued measure 
(POVM) of on-off detector is described as 
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FIG. 2: Probability distribution of the single-mode non- 
Gaussian state (sòlid line, A = 0.4, T = 0.9) and the squeezed 
vacuum state (dotted line, A = 0.4) for the phase parameter 
ip — with ideal setup. 



FIG. 1: Measurement-induced non-Gaussian operation on 
the single-mode squeezed vacuum state. BS, PD, M, and HD 
are beam splitter, photon detector, mirror, and homodyne 
detection, respectively. 



Simultaneous on events in both modes C and D project 
the state over mode A and B into the state 
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is the success probability oí this on event selection, and 
A = tanh r. Finally, remaining modes A and B are recom- 
bined with another beam splitter of the transmittancc 
t(= 0.5), being transformed as 
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The output state of mode B is the vacuum state due to 
the interference. 

The success probability and the quality of the output 
non-Gaussian state can be controlled by the transmit- 
tance of the tapping beam splitter T. There is a trade- 
off between these two. The non-Gaussian property ap- 
pears more strikingly for higher transmittance but with 
sacrifice of lower success probability. Hereafter, we set 
T = 0.9. 

The output non-Gaussian state in mode A is measured 
by the homodyne detection. The measurement basis is 



given by the the quadrature eigenstate 
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The probability distribution is given by 
Phd(x v ; X) 



where 

1 



Pn(x v ;X) - P w {x v ;X) 

-Poi{x v ; X) + P 00 {x v ; X),{13) 
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and 7n = R = 1 - T, 710 = 701 = R/2, and 700 = 0. 
Thus the probability distribution consists of the four 
Gaussian terms. Fig. shows the probability distri- 
bution of non-Gaussian state (A = 0.4) along x v= o = x 
axis (the sòlid line) and the one of the input squeezed 
vacuum state (the dotted line). The probability distri- 
bution of the output non-Gaussian state consists of the 
single main peak with two small side lobes. The prob- 
ability distribution dominated by the single main peak 
may still be characterized by the variance 
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FIG. 3: Variance of the single-mode non-Gaussian state (sòlid 
line, T = 0.9) and the squeezed vacuum state (dotted line) 
with ideal setup. 
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and Jij's have been given above. In Fig. |3 the variances 
of the output non-Gaussian state and of the original input 
squeezed state are compared. As seen, in the range of 
A < 0.47, the variance of the output non-Gaussian state 
is lower than that of original squeezed state, which mcans 
that the squeezing degree is effectively enhanced. 

Our interferomctric non-Gaussian operation scheme is 
actually found to be equivalent to the Schròdinger cat- 
like state generation proposed in Ref. 01 (Fig- EJ- This 
can directly be observed by the following equation 
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FIG. 5: Wigner function of the single-mode non-Gaussian 
state (A = 0.4, T = 0.9) with ideal setup. 



by noting that input modes B, C, and D are all the vac- 
uum state. 

The Wigner function 
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of the non-Gaussian state is obtained by 
Wng(x,p;\) 
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This is illustrated for two A's in Figs. [SJ (A = 0.4) and|Ü] 
(A = 0.8). For smaller A, the Wigner function is closer 
to that of the squeezed vacuum state. For larger A, on 
the other hand, the property of cat-like state becomes 
more remarkable. These figures will be compared later 
with the case including imperfections. 
Interestingly, the average photon number 
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increases after conditional non-Gaussian operation (Fig. 
0), where N = à^à is the photon number operator, and 
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FIG. 6: Wigner function of the single-mode non-Gaussiari 
state (A = 0.8, T = 0.9) with ideal setup. 
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FIG. 7: Average photon number of the single-mode non- 
Gaussian state (sòlid line, T = 0.9) and the squeezed vacuum 
state (dotted line) with ideal setup. 



III. NON-GAUSSIAN OPERATION ON 
TWO-MODE SQUEEZED STATE 

Now we turn to the two-mode case (Fig. |HJ in the ideal 
situation, where the other squeezed vacuum state with 
the reduced uncertainty of the p (= )-quadrature 
is input into mode B instead of the vacuum state. The 
two-mode squeezed vacuum state 
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is generated after the first 50:50 beam splitter, where 
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(24) 



is the two-mode squeezing operator. The operations after 
this is similar to the single-mode case. The state condi- 
tioned by the simultaneous on results for modes C and 
D is 
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FIG. 8: Measurement-induced non-Gaussian operation on 
the two-mode squeezed vacuum state. BS, PD, M, and HD 
are beam splitter, photon detector, mirror, and homodyne 
detection, respectively. 



where 
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with the high transemittance T — cos 2 9, and 
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is the success probability of the on event selection. 

The two modes of the conditional non-Gaussian en- 
tangled state are combined via the second 50:50 beam 
splitter, and are then measured by the two homodyne 
detectors, each of which measures the two orthogonal 
quadratures simultaneously. This chain of operations 
corresponds to the CV Bell measurement represented by 



|n(x,p)) fcí 



dyé^\x + y) k \y)i. (28) 



Therefore, the probability distribution of homodyne de- 
tection in phase space is 



n(x,p)|^.|n(a;,p)) 

P£\x,p;\)+P$(x,p;\), (29) 
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FIG. 9: Probability distribution of the two-mode non- 
Gaussian state (sòlid line, A = 0.4, T = 0.9) and the squeezed 
vacuum state (dotted line, A = 0.4) for the phase parameter 
tp — on the mode A with ideal setup. Probability distribu- 
tion for tp — 7r/2 on the mode B gives the same result. 
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and 7Í 2) = R, jg* = 0. We can see that Eq. 
is also the combination of four Gaussian terms, which 

is totally non-Gaussian distribution. Fig. [ÜJshows the 

(2) 

probability distribution P^(x,p; A) along x axis for A = 
0.4, after integrating out the variable p. In this two mode 
case, the two small side lobes cannot be seen, and the 
distribution is close to Gaussian. This distribution may 
be characterized by its variance, 
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This is shown in Fig. Eïlby the sòlid line. The dotted line 
corresponds to the case of the two-mode squeezed vac- 
uum state. In A < 0.67, the variance of two-mode non- 
Gaussian state is lower than that of two-mode squeezed 
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FIG. 10: Variance of the two-mode non-Gaussian state (sòlid 
line, T = 0.9) and the squeezed vacuum state (dotted line) 
with ideal setup. 



vacuum state. Thus the squeezing degree is effectively en- 
hanced in the two-mode case as well. The mcasurement 
on CV Bell basis <|28|) is one of the correlation between 
two modes, therefore we can consider the variance <|3l·ll as 
a measure of quantum correlation of bipartite entangled 
pair. 



IV. APPLICATION TO DENSE CODING 
SCHEME 



The entanglement measures for CV systems have been 
clarified for Gaussian states so far. But the measures 
for mixed non-Gaussian states are not clear yet. In this 
regards, it might be sensible to adopt some operational 
measures connected directly to CV protocols. In this 
scction, we evaluate the property of the conditional non- 
Gaussian state in terms of the mutual information for 
the dense coding scheme proposed in 0, [2^ ■ The result 
here is added to the ones based on the other kinds of 
operational measures. such as the improvement of the 
teleportation fidelity [lg and the nonlocality due to the 
violation of Bell type inequality 0, 0, 0, [2(j • 

The variance reduction of the non-Gaussian states 
shown in the previous section implies the improvement 
of the signal-to-noise ratio, and hence the increase of the 
transmissible information. This can be measured by mu- 
tual information, which is calculated by evaluating the 
channcl matrix between the input and output signals. 
This measure is known to be a stringent figure of mèrit 
for the system in the sense that the gain usually vanishcs 
even with relatively small imperfections. Therefore it 
would be a nice criterion to see whether the information 
gain is observed or not by the non-Gaussian operation 
compared with the original two-mode squeezed vacuum 
state. 

In Fig. ^2 the setup of dense coding scheme is de- 
picted. Alice encodes classical message (x s ,p s ) on one of 
the beam of the two-mode non-Gaussian entangled state 
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by the displacement operation, 



U A (x s ,p s ) = e-**«"D 



X s + ÍPs 

V2 



(33) 



where D{a) = exp(aa^ — a*à). Now, we consider a sim- 
ple model in which signals x s and p s are restricted to 
±\/2Q! with equal prior probabilities, (quadrature phase 
shift keying) and thus a set of 2-bit information {üqq = 
(x s = V%a,p s = V2a), a Q1 = (V2a, — s/2a), a xo = 
(—V^a,V2a), a\\ = (— V2a, — V2a)} is encoded. (The 
continuous encoding onto the non-Gaussian entangled 
state makes the calculation difficult, so we have here sim- 
plified to the 4 valued discrete encoding. ) 

Bob attempts to decode classical message from Alice 
(x s ,p s ) by the CV Bell measurement described in Eq. 
(|2*S|) . The distribution of the detection probability is 
given by 



{n(x,p)\Ü A (x s ,p s )p^' G Ü A (x s ,p s )\IL(x,p)) 
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where we use the relation 
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Then Bob identifies the 2-bit classical messages according 
to the decision rule {&oo — {x > 0,p > 0), &oi = (x > 
0,p< 0), b 10 = (x<0 7 p> 0), b n = {x<0,p< 0)}. 

Now the 4-by-4 channel matrix is given by 
{P(a mn \bki)} where the elements for (m = k,n = l) con- 
tribute to the success probability and otherwise, the error 



probability. For example, 
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is the error function. Other components can be obtained 
similarly, 
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x (l + (-l) m - fe crf[fly(A)a]) 

x (1 + (-l) n -'erf [%(A)a]) . (39) 

The mutual information is calculated by the above chan- 
nel matrix as 

J(A;B) 

= ^ P (aki)P(b mn \aki) 



i 2 ^ 



x log 2 
2 [bit] 



P{b mn \a k i) 



^P(a M )P(b \ctki) 

k,l 

as A — > 1, 



(40) 



where P{aki) = 1/4 is the prior probability. 

From the viewpoint of communication efficiency, we 
should optimize the mutual information under the power 
constraint condition that the total energy of initial 
squeezing power and encoding displacement power is con- 
stant. We are, however, interested in how to quantify 
the improvement brought by the entanglement enhance- 
ment due to the non-Gaussian operation for a given in- 
put. What should be fixed is then the degree of entan- 
glement of the input state, equivalently the squeezing 
degree A for the input. We intend to compare the mu- 
tual informations between the cases with and without the 
non-Gaussian operation. 
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Squeezed vacuum state 
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X 

FIG. 12: Mutual information via the dense coding channel 
(a = 1.5) of the non-Gaussian state (sòlid line, T = 0.9) and 
the squeezed vacuum state (dotted line) with ideal setup. 
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cc=0.7 



Non-Gaussian state 
Squeezed vacuum state 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

X 

FIG. 13: Mutual information via the dense coding channel 
(a = 0.7) of the non-Gaussian state (sòlid line, T = 0.9) and 
the squeezed vacuum state (dotted line) with ideal setup. 



In Figs. ll2land l*Hfl the mutual information of a = 1.5 
and a = 0.7 are indicated respectively, where the orig- 
inal squeezed vacuum state cases are also indicated by 
the dotted line. We can see triat, in the range of 
A < 0.38 and A < 0.65, in Fig. [Ï3and[ï31 respectively, 
the mutual information of two-mode non-Gaussian state 
is larger than the original two-mode squeezed vacuum 
state case. This can naturally be regarded as the gain as- 
sistcd by the enhanced entanglement of the non-Gaussian 
cntangled state. Via the probabilistic photon subtraction 
operation, similar to the single-mode case, the average 
photon number of output two-mode non-Gaussian state 
increases compared with input squeezed vacuum state, 
which brings effective enhancement of entanglement. 



V. POVM OF PRACTICAL PHOTON 
DETECTOR MODEL 

From this section, we consider possible imperfections, 
including the linear loss in optical paths, finite quantum 
efhciency and nonzero dark counts of the on-off detector. 
The imperfect photon detector may be modeled by the 
following POVM element [Ü, 



n(m;„)= £ ( Jn )r> m 0--V) n - 



\n)(n\ 



(41) 



where n photons are converted to m < n with the imper- 
fect quantum efhciency rj, and 



(n — m)!m! 



(42) 



is the binomial coeficient. This element is considered 
as the sequence of two devices; beam splitter of trans- 
mittance r\ (linear loss) and perfect n photon counting 
detector, where n — m photons are removed as the linear 
loss. As a whole, m photons are detected. 

Further, assume that the net count of photon detector 
is N(> m), where N — m photons are due to dark count, 
following Poisson distribution. With average dark count 
v, the POVM of detector counting TV photons is 



n(/v,77,i/) 

TV 

= £< 

m=0 



,N-m 



(N — m)\ 



E 



r, m {l-rf) n - m \n){n\, 



which satisfies completeness relation, 

00 

£n(i\r, »?,!/) = i. 



(43) 



(44) 



N=0 



where 



With Eqs. and gU, we obtain a set of POVM's of 
practical on-off type photon detector, 

{n^^n^M}, (45) 
n(°%,i,) = n(o,77,^) 

00 

= e -^(l-77)»(n|, (46) 

00 

n (on W) = 53 n(jv, »/,!/) 

N=l 

= i-fi^(n,y). (47) 



VI. ANALYSIS WITH PRACTICAL 
PARAMETERS IN SINGLE-MODE SCHEME 

Let us first go with the single-mode non-Gaussian oper- 
ation (Fig. The quantum efuciency of the homodyne 
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detector may be considcrcd as approximately the unity, 
which is the case of the time limited signals prepared by 
chopping the contirmous wave field generated from the 
cavity enhanced optical parametric process. We assume 
that the linear loss occurs between the first and second 
beam splitters, which may be attributed to the loss on 
beam splitters, mirrors, nonlinear crystals. Such a de- 
cohered squeezed vacuum state is transformed into the 
non-Gaussian state by the imperfect on-off type photon 
detectors modeled above. 

The linear loss in optical paths can be described with 
a beam splittcr of transmittance Tl, 

IVOABCDL1L2 

= VÍd(0)Vac (0)^(0 VAL t (0 

xVab (I) |r) A |0) B CDL 1 L 2 (48) 

where L\ and L 2 are assigned for loss channels and 



tan£ = 



l-T L 



(49) 



Modes of Li and L 2 are traccd out, and mode C (D) 
is tapped from mode A (B) by a beam splitter of high 
transmittance T respectively. Therefore input state is 

Qin = Tr LlL2 [K n )(ABCDL 1 L 2 )(V'ínl] ■ (50) 

Then, simultancous on evcnts on modes C and D are 



selected conditionally, 



Tr 



CD 



nW(^)®n7(^) 



(on) . 



£>out = 

where 

Vdct = TrABCD 

= l-2e~" 



Vdc 



éi» ® ® ïi^ 1/) 



(51) 



(on) , 



+e 



-2v 



1 - À 2 

y 1 - \ 2 (T L T + R L + ^T L R) 2 

1 1-A 2 
1 - \ 2 {T L T + R L + (1- r,)T L R} 2 

(52) 



is the success probability of event selection (Rl = 1 — 
Tl). Then we recombinc modes A and B with another 
balanced beam splitter and mode A is measured with 
homodyne detection (mode B is vacuum state), 

= 'Pu(x lp ;X) — e~ v Vw(x íp \ A) 

—e~ l "P m {x ip -.\) + e~ 2v V 00 (x lp ; A), 



(53) 



where 



Í?NG 

{xtp , A) 



Mi 

i 



Qout 



(54) 



1-A 2 



V5FPdet V (1 - AT L T) 2 - A 2 (i? L + V-) 2 + 4AT L Tsin 2 <p 



x exp 



l-X 2 (T L T + R L +i l3 ) 2 



(1 - AT L T) 2 - A 2 (iï L + V.)2 + 4AT L Tsin 2 <p v 



(55) 



and y n = T L R, i w = j' 01 = (2 - r 1 )T L R/2, and where 

7oo = (1 — v)TlR- With above results, we can obtain 

the variance of non-Gaussian state along x axis (ip = 0), V (A) 



1-A 2 



27>dct V 1 - A 2 (T L T + R L + 7-j) 2 
(l-AT L T) 2 -A 2 (i? L+7í ' ? ) 2 



x ■ 



1-X 2 (T L T + R L +^ 



(57) 



We assume the total linear loss 25% (Tl = 0.75), and 

V(A) 

the on-off type detector of quantum emeieney r\ = 0.6 
and the dark count rate 10000 [counts/sec]. At present, 



/°° f r°° 1 and the dark count rate 10000 [counts/sec]. At present, 

dx x Vut>(x\ A) — < / dx xVhd(x; A) > lower dark count rate ~ 100 [counts/sec] is achievable in 

"°° _ v _^ ~°° laboratory, however, we consider the rather large value 

= Vn(A) - e "Vio(A) - e "Voi(A) + e "Voo(A), to see the effect more clearly. When the gating time of 

(56) photon detector is 10~ 7 [sec], the net dark count in a 
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Single-mode 

l 2 -(Practical) 
1 - 



Non-Gaussian state — 
Squeezed vacuum state 

1=0.4, r= 



W(x,p;X) 



0.9 




FIG. 14: Probability distribution of the single-mode non- 
Gaussian state (sòlid line, A = 0.4, T = 0.9) and the squeezed 
vacuum state (dotted line, A = 0.4) for the phase parameter 
ip — with practical setup (Tl = 0.75, r\ = 0.6, v = 10~ 3 ). 



Non-Gaussian state 
Squeezed vacuum state 



T=0.9 



5 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

A 

FIG. 15: Variance of the single-mode non-Gaussian state 
(sòlid line, T — 0.9) and the squeezed vacuum state (dotted 
line) with practical setup (Tl = 0.75, r) = 0.6, v = 10~ 3 ). 



single event is v = 10 3 , which is large enough to affect 
the output state. 

Figure El shows the probability distribution of non- 
Gaussian state (A = 0.4) along x axis (sòlid line) and the 
one of the input squeezed vacuum state (dotted line), 
where the two side lobes seen in Fig. [21 disappear due 
to imperfections. In Fig. El the variances of the output 
non-Gaussian state and of the original input squeezed 
state are compared. The reduction of the variance due 
to the non-Gaussian operation is seen for A < 0.4. At 
this cross point, the variance shows 2.5 dB (in terms 
of -101og 10 [V(A)/V(0)] [dB]) below the shot noise level. 
The range of the variance suppression becomes narrower 
than the ideal case where A < 0.47 (Fig. [3J. This is 
mainly due to the dark counts. The linear loss spoils the 
degree of squeezing, as the increase of the overall variance 
level. The imperfect quantum efhciency reduces success 
probability Vdct- 

The Wigner function of the non-Gaussian state in the 




FIG. 16: Wigner function of the single-mode non-Gaussian 
state (A = 0.4, T = 0.9) with practical setup (T L = 0.75, 
T) = 0.6, v = 10~ 3 ). 



W(x,p;X) 




FIG. 17: Wigner function of the single-mode non-Gaussian 
state (A = 0.8, T = 0.9) with practical setup (T L = 0.75, 
ti = 0.6, v = 10~ 3 ). 



imperfect setup is calculated as 



Wng(^,p; A) 
= A) - e-"W w {x,p; X) 



-e- v Wox(x,p; A) 



'Wooix^p; A),(58) 



where 



i 



l-A 2 



T^Vdc 



x exp 



x exp 



l-V{T L T-R L -i l0 Y 
l-X 2 {T L T + R L +y tJ ) 2 



(1 - \T L Tf 



{l-\T L Tf-\ 2 {R L +j t] f 2 
l-A 2 (T L T-iï L - 7l ' J ) 2 P 



(59) 



The Wigner functions for A = 0.4 and A = 0.8 are illus- 
trated in Figs. lluland lT7l respectively. The former is still 
similar to the squeezed state, while the latter is close to 
decohered plus-cat state, compared with Figs. [SI and 
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VII. ANALYSIS WITH PRACTICAL where 
PARAMETERS IN TWO-MODE SCHEME 



We turn to the analysis of two-mode scheme with prac- 
tical parameters (Fig. |8J. Similar to the single-mode 
case, we assume that the linear loss occurs between the 
first and second beam splitters, and model the effect by 
inserting a beam splitter of transmittance Tl into each 
arm. The resulting state is then 



where 



~(2) 



//(2)\ 



I /'( 2 )\ 

m n i (ABCDLiL 2 ) 



'(2) 



(60) 



IV'in ' /AB< 'DL,L_. 

= Vt 



x |r) A | - r·) B |0)cDL 1 L 2 · 



(61) 



The conditional state selected by simultaneous on events 
on modes C and D is then, 



V. 



(2) 
del 



1 - 2e~ 



1 - A 2 



-2v 



l-A 2 {l-r/T L (l-T)} 
1 - A 2 



1-A 2 {1-77T L (1-T)} 2 



is the success probability of this event selection. 
probability distribution of the Bell measurement is 



(63) 



The 



3 (2) 
HD 



(x,p; A) 



U(x,p)\g^\TÍ(x,p)) 



v[f(x,p; \)~e-»Viï(x,p; A) 

-e-"V$ (*> P5 A ) + e ' 2u Ho (*> P! A) , (64) 



— i"n(2) / 



.(2) 



Tr C D 



n 



(on) 



n 



(on) 



V. 



(2) 
det 



(62) 



where 



-P^(x,p;X) 



1-A 2 



2nV£\ (1 - AT L T) 2 - A 2 (i? L + 7 f } )(í?l + 1 3 



'(2h 



x exp 



1 - A 2 (T L T + R L + 7 ; (2) )(T L T + R L + 7 ' (2) ) ' 

-(x z +p z ) 



2{(1- AT L T) 2 - A 2 (i? L + 7 f >)(i? L +7D) 



'( 2 h 



(65) 



and 7Í l } = T L R, -/ (Z> = (1 - r/)T L R. With Eqs. © and we can calculate the variance of the output state at 
each port, 



V< 2 >(A) = vg>(A) - e-"VÍ?(A) - e^V ( ?(A) + e^V^A), 



(66) 



where 



'dot l 1 



1 (1 - A 2 ){(1 - AT L T) 2 - A 2 (i? L + ^ 2) ){R L + tD) 



A 2 (T L T + R L + 7 - (2) )(T L T + R L + 7 '^)} 2 



'(2)n 



(67) 



Fig. H8lshows the probability distribution of the output 

( 2) 

non-Gaussian state Vl l f ) (x,p; A) along x axis for A = 0.4 
(sòlid line) , after the integrating out the variable p. We 
can see that the peak is sharper than the input squeezed 
vacuum state of A = 0.4 (dotted line). The variance of 
the output non-Gaussian state V(A) is shown in Fig. ^3 
Its variance becomes lower in the range of A < 0.63 than 
that of the input squeezed vacuum state (dotted line), 
which corresponds to the range up to 3.8 [dB] in the 
relative decibel scale. 



Now, let us consider the dense coding scheme with the 
two-mode non-Gaussian state (|60f) in accordance with 
Sec. IIVI We assume that the decoherence of the en- 
coded signals in the way from Alice to Bob is negligible. 
Actually we are interested in quantifying the effect of 
non-Gaussian state in terms of the mutual information 
rather than in applying the scheme to practical Commu- 
nications. The procedure goes in a similar way to the 
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Two-mode 
0.8 (Practical) 



Non-Gaussian state 
Squeezed vacuum state 



A=0.4, 7=0.9 




ideal case (Fig. HT|l . 



n(a; , p) | U A (x s , Ps ) £?ng U A ( x »> Ps ) | n(z , p) ) 



p££(a;-x„p-p > ;A).(68) 



FIG. 18: Probability distribution of the two-mode non- 
Gaussian state (sòlid line, X = 0.4, T = 0.9) and the squeezed 
vacuum state (dotted line, A = 0.4) for the phase parame- 
ter íp — on the mode A with practical setup (Tl = 0.75, 
77 = 0.6, v — 10~ s ). Probability distribution for ip — ty/2 on 
the mode B gives the same result. 



With Eq. I|68|) , we can obtain the channel matrix similar 
to Eq. (|36|l and the mutual information can be obtained, 



Non-Gaussian state 
Squeezed vacuum state 



7=0.9 - 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

X 

FIG. 19: Variance of the two-mode non-Gaussian state (sòlid 
line, T = 0.9) and the squeezed vacuum state (dotted line) 
with practical setup (Tl = 0.75, 77 = 0.6, v = 10 -3 ). 



J(A;B) 



^ P{aki)V(b mn \a k i) 



k,l,m,n 



x log 2 



T(b 

\am) 



\au) 



(69) 



where 



V(b m n\aki) 



E 



1 - à 2 



^0 1 - A 2 (T L T + R L + ir)(T L T + R L +1 3 



' (2) )(T L r + i? L + 7 f ) ) 
x (1 + (-l) m - fc erf (X)a\ ) (l + (-l)"-'erf (X)a] ) 



(70) 



and 



ÜUX) 



1 - X 2 (T L T + R L + l[ {2) ){T L T + R l + j^>) 
(1 - XT L TY - X>(R L + 1 ' l (2) ){R L +^ 2) ) 



\ 



(71) 



In Figs. I2UI and ÜH1 the mutual information of a = 1.5 
and a = 0.7 arc indicated, respectively. The mutual 
information of the non-Gaussian state increases over that 
of the original squeezed vacuum state, in the range of 
A < 0.47 for a = 1.5 and in A < 0.65 for a = 0.7. 
This gain may be attributed to the effective increase of 
the entanglement, and will be seen even under practical 
situation. 
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FIG. 20: Mutual information via the dense coding chan- 
nel (a = 1.5) of non-Gaussian state (sòlid line, T = 0.9) 
and squeezed vacuum state (dotted line) with practical setup 
(T L = 0.75, T) = 0.6, v = 10~ 3 ). 



2 - 



1.8 




FIG. 21: Mutual information via the dense coding chan- 
nel (a — 0.7) of non-Gaussian state (sòlid line, T = 0.9) 
and squeezed vacuum state (dotted line) with practical setup 
(T L = 0.75, 77 = 0.6, v = 10" 3 ). 

VIII. DISCUSSION AND CONCLUSION 

In this paper, we have studied the non-Gaussian oper- 
ations induced by the measurement with the on-off de- 
tectors on the single- and two-mode Gaussian squeezed 
vacuum states. Our scheme is the Mach-Zehnder inter- 
fcrometer, where the non-Gaussian state measured at the 
two output ports with respect to the single quadrature 
in each port. This setup was originally motivated by 
the naive intuition that the useful non-Gaussian states 
induced by the measurement should exhibit smaller vari- 
ances in appropriate quadratures than the ones of the 
originals, irrespective of the mixedness of the state and 
the amount of the deviation from the Gaussian state. 
In both the single- and two-mode cases, the homodyne 
probability distributions at the two output ports show 
the single main peak which is very close to the Gaussian. 
The effect of the non-Gaussian operations based on the 



on-off detector and the usefulness of the resulting states 
may be characterized simply by the reduction of the vari- 
ances. In the two-mode case, the two variances measured 
at the two output ports can still be regarded as a measure 
of the quantum correlation in the non-Gaussian bipartite 
system induced by the on-off detector. This sort of evo- 
lution is of course not rigorous mathematically. More 
satisfactory theories need to be developed. 

Our analysis includes possible practical imperfections, 
such as the linear loss in optical paths, finite quantum 
efhciency, and nonzero dark counts of the on-off detector. 
We have seen that although they degrade the output non- 
Gaussian state, the gains in terms of the variances and 
the mutual information of the dense coding scheme still 
remains. One of the important effects that has not been 
involved in this paper is the mode mismatch between 
the field measured by the homodyne detector and the 
field of trigger photons measured by the on-off detector. 
Ideally, the field mode of trigger photons must be in a 
mode which overlaps perfectly the local oscillator mode 
in the homodyne detector, with respect to the spatial, 
temporal, and frequency domains. Otherwise one cannot 
select the matched mode, and is led to the degradation of 
the gains. This aspect was studied in [25|. An alternative 
analysis based on the present model will be presented 
elsewhere. 

Then, as an operational measure of non-Gaussian oper- 
ation, we have studied the mutual information in quan- 
tum dense coding scheme. The dense coding is one of 
the entanglement-assisted schemes, and the improvement 
of the mutual information is considered as the gain as- 
sisted by the enhanced entanglement via non-Gaussian 
operation. The mutual information in the dense coding 
scheme has some advantageous features. Firstly, it is 
a well-established measure in information theory, which 
is obtained by specifying the channel matrix, and has a 
clear operational meaning in a múltiple use of the chan- 
nel. Secondly, it can be a rigorous measure for the ob- 
tained gain. Actually, the gain in the mutual information 
usually gets lost even by a small amount of imperfections. 
So the observed gain will ensure that the system really 
works. Finally, to evaluate the mutual information is 
also suitable for practical use. In fact, the corresponding 
experimental setup is relatively easy in laboratory com- 
pared with, for example, quantum teleportation. This 
measure would be the first quantity to evaluate when the 
non-Gaussian state is experimentally available. 

On the other hand, more strict theories for measure 
of mixed entangled state has to be developed, which is 
earnestly desired both theoretically and experimentally. 
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